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Abstract 

As scientists we like to think that modern societies and their members base their views, opinions and 
behaviour on scientific facts. This is not necessarily the case, even though we are all (over-) exposed 
to information flow through various channels of media, i.e. newspapers, television, radio, internet, and 
web. It is thought that this is mainly due to the conflicting information on the mass media and to 
the individual attitude (formed by cultural, educational and environmental factors), that is, one external 
factor and another personal factor. In this paper we will investigate the dynamical development of opinion 
in a small population of agents by means of a computational model of opinion formation in a co-evolving 
network of socially linked agents. The personal and external factors are taken into account by assigning 
an individual attitude parameter to each agent, and by subjecting all to an external but homogeneous 
field to simulate the effect of the media. We then adjust the field strength in the model by using actual 
data on scientific perception surveys carried out in two different populations, which allow us to compare 
two different societies. We interpret the model findings with the aid of simple mean field calculations. 
Our results suggest that scientifically sound concepts are more difficult to acquire than concepts not 
validated by science, since opposing individuals organize themselves in close communities that prevent 
opinion consensus. 



Introduction 

Human societies are examples of dynamical complex systems where the complexity stems, on one hand, 
from the social interactions between individuals of the society, and on the other from societies' organizing 
structure, making them challenging to study empirically. In these systems the network approach and 
mathematical modelling have turned out to be crucial in gaining understanding of their structure and 
dynamical behaviour [lj|3] . This has been achieved by using such techniques as the analysis of large data 
sets, collected in questionnaires or obtained by automatically recorded digital footprints of our actions. 
One of the interesting phenomena in a society - describable as social networks of individuals and links 
between them - is the process of opinion formation among people, where it is assumed that people focus 
on a single question answerable simply by yes or no. Likewise one can consider a rather similar situation 
of spreading scientific information among people, such that during the process some people adopt it as 
truth and some others reject it. To address these issues there has recently been increasing interest in 
modelling opinion formation, with many examples in the literature |4|-(9]. Here we extend this approach 



to study the response of people to scientific facts in co-evolving social networks 10 ■ 14 and under the 



influence of mass media 15 ■ 17 



In a recent paper [18| we have described a new model for opinion formation, based on co-evolution 
dynamics with a single continuously-changing state variable (measuring the degree of agreement with a 
single question posed) and with link rewiring, i.e. changing the structure of the network. These two 
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mechanisms take place with different time scales such that the former stands for a fast transaction time 
scale and the latter for a slow generation time scale. The rewiring mechanism is based on simple rules 



adopting the social network dynamics proposed by Kossinets and Watts 19 , that is, a drive for the 
closure of network triangles and similarity of individuals' opinions. Our model was designed to capture 
some of the most essential features of opinion formation, namely, dyadic interactions corresponding to 
one-to-one discussions, an attitude parameter describing the personal reaction of an individual to the 
overall opinion, and an external influence on the whole network of individuals mimicking the effect of 
media. In our previous studies 18 20 we ignored the role of an external field to simplify the treatment 
of the model. In reality, however, it is difhcult to envisage a situation in which the external field is 
zero, since there is always information available to all individuals that can affect their decision-making 
processes. 

A good example of this situation is how much scientific and technological knowledge determines the 
public engagement with polemic issues that affect society. Although we live in societies deeply intertwined 
with modern science and various emerging technologies, most of the world's population is scientifically 



illiterate, as can be inferred from examining children of school age 21 . In contrast one would expect 



a literate citizen to be able to evaluate the quality of scientific information, take rational positions on 
polemic decisions, and make a distinction between science and pseudo-science. A very important goal of 
science literacy and its public understanding is precisely to erase scientific mistakes and superstitions ^22, . 
In order to understand the public perception of science and technology different models have been 



built, but there is still some controversy between the two main approaches 23 . First, the science literacy 



model assumes that knowledge increases the public acceptance of scientific research and allows scientists 
to work with more freedom. As a consequence a knowledgeable public would reach consensus with experts 
quite easily. The second is a "contextualist" approach to the public understanding of science, where social 
values and attitudes of the public are also very important. This is the so-called ethnographic case study 
method, where alternative forms of knowledge interact with the actions of experts. Quantitative surveys 
are often related to the science literacy model and generalized attitudes, while the ethnographic approach 
emphasizes local knowledge. This conflict has led to the dismissal of either surveys or ecological data, 
although the two models can complement each other. In 23 a survey-based approach was proposed to 
integrate multiple influences and bridge this gap. 

In a recent study 24 a comparison between roughly 200 surveys from 40 countries was made, find- 
ing a correlation between knowledge of scientific facts and a favourable attitude of the public towards 
science. However, it also turned out that knowledge itself seems to have just a small influence on the 
opinion. Instead moral values, religious beliefs, trust and cultural trends of communities were found to 
be significant factors for building the final attitude of an individual. This study concluded that more 
research needs to be done on the relationship between knowledge and attitude towards science. 

It is commonly considered desirable that the development of science in a country should be linked to 
a policy of science communication, leading to an increase of the public understanding of science. The 
term "science communication" is used here as a technical term to describe the multidisciplinary field of 
science dissemination. In a recent book on the matter 25 the editors report several stories of success 
in science communication, covering a wide range of topics in which this link is associated with success. 
They also point out that the scientific community is increasingly aware of the importance of such a 
link. However, the activities related to communicating science to a broad audience are not based on 



quantitative models, and the main source of data is still heavily based on surveys 26 . In present day 
surveys very little attention is put into extracting dynamical changes in the general knowledge of the 
public, which are extremely important to testing the success of different policies. 

In general, there is a great interest in how people react to new technologies, how concerned citizens are 
about their risks, and how aware they are of their benefits. For instance, in the case of nanotechnology a 
recent survey in USA [27] found the opinion of the public to be quite favourable. This study shows that 
if there is an effort to give simple and accurate information the public will react in a positive way, even 
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enabling them to associate nanotechnology with other fields such as medicine, machinery or computing, 
and thus tending to have positive opinions and a perception of low risk. 

As a summary of the research performed towards the public understanding of science, we list the main 
factors known to influence social opinion: 



Knowledge. In 24 a weak correlation between knowledge of scientific facts and positive attitude 
towards science was found, but also that it can become negative in case of certain technologies, such 
as in the human embryo research. These results show that people who are more scientifically literate 
have more positive attitudes towards science in general, but can be against a specific application or 
scientific research. 

Science Communication. In relation to this factor the quality and relevance of science communica- 
tion information could move opinion from negative to positive, as in the case of nanotechnology [27' . 



• Culture. When different countries are compared social, political, and cultural disparities between 
nations affect the results in fundamental ways. Local ideology is yet another factor. 

All these factors combined affect the personal opinion. There is also a relation with regional socio- 
economic conditions: in the developing countries science is often idealized, while in post-industrial soci- 



eties there is not such a positive cultural stereotype 24 



In addition to the science perception models discussed above, a remark must be made on cognitive 



models. Many of these 23 27 are based on the cognitive miser model, where citizens are supposed to 
collect only as much information about the topic as they think is necessary to reach an opinion ^28^ . Most 
of the citizens avoid deep information but rely more on their ideology together with information from 
mass media. The news, reported in different media, is the main component for the opinion formation 
process. The cognitive miser model of social thinking reflects the importance of motivation and emotions 
and the reluctance to do much extra thinking. Besides this superficial associative reasoning there is a 
deeper cognitive process, which can be considered true reasoning. The mind has two different modes of 
processing, one automatic and another conscious, which in a way work together. Some situations require 
the individual to overcome unprecedented obstacles, as happens in scientific research, and the conscious 
mind takes the challenge. This explains the attitude of experts towards certain problems. As for the two 



cognitive processing modes 29 , the average citizen is usually in the cognitive miser state most of the 
time. 

Considering all these factors and cognitive models, the main motivation of the present work is to 
use the knowledge gathered by scientists in different fields, including social scientists, to construct a 
mathematical model that could lead to the quantification of the process of public understanding of 
science. With this information in mind we have extended our opinion- formation/co-evolution model 
to simulate the effects of external information on perception spreading among individuals in a social 
network. The key ingredients are the following: Each citizen has a personal opinion that is related not 
only to his or her scientific education but also to the cultural background, and that is likely to change 
through personal one-to-one interactions. There is an external influence coming from media and science 
communicators represented by a field. Finally, as the long-range order and interactions are also important 
and the perception of the general information is considered personal, they are represented by an attitude 
parameter that corresponds to the reaction of an individual to what friends, relatives, and colleagues 
think. 

This paper is organized such that first we describe the model we have used to simulate the effect of 
external bits of information on the opinion formation in a social network of individuals. Then we study 
the model's behaviour via numerical simulations and analytical approximations. After that we adjust the 
field strength by using data from two survey studies realized as questionnaires. Finally we discuss our 
results and present concluding remarks. 
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Results 

Description of the model 

As our model to simulate the effect of external information coming from media on the formation of 
perception of an issue in a social network of interacting individuals, we take the co-evolving opinion 
formation model first studied in [18) . This model is used since we consider it to be quite generic in 
nature, and straightforwardly extendable to include the influence of media by introducing an external 
field that affects all individuals of the society. Thus the equation for the state variable of each agent i 
can be written as, 

"dt^ = (1) 

j 

where the operator O stands for the rewiring operations performed on the adjacency matrix Aij of the 
social network describing the pairwise interactions between agents i and j. These changes occur after g 
opinion transactions have taken place, defining a slow time scale T = gdt for the network evolution. Each 
agent is characterized by a state variable Xi representing the instantaneous inclination of the personal 
perception or opinion on a posed question or issue, and it is bounded between -1 (total disagreement) 
and 1 (total agreement). The fast dynamics of each state variable is described as follows. 



~dt 



(^ifl{{Xj}l) + Xifs{{Xj)s) + hi, (2) 



where the first term on the right-hand side represents the reaction of agent i to the overall opinion of 
all distant individuals {xj}i, modulated by the agent's own attitude (a^) towards the overall or public 
opinion Hence we can write, 

where mi{i) means the set of nodes that are £ steps away from node i, and ^max is the number of steps 
needed to reach its most distant neighbours. Note that this quantity is different for each agent, since 
it depends on the local topology of the network. The second term on the right-hand side of Eq. ^ 
represents discussions between pairs of linked agents, such that, 

/s=sgn(xj) ^ Xj. (4) 

The resulting absolute value \xi \ in Eq. ([2| assures that agents with the same sign of opinion reinforce their 
convictions after a discussion, while agents with opposite sign become less convinced of their position. 
These are considered as reasonable assumptions based on simple social interactions. The last term hi on 
the right-hand side of Eq. ([2]) is an external field representing the personal bias towards either opinion 
(-1 or -l-l), due to e.g. mass media (newspapers, TV, radio), and will be discussed further at the end of 
this section. 

Note that \xi\ in Eq. ^ could eventually become larger than one, which we choose to interpret as a 
state of total conviction. Thus the x variables of totally convinced agents are set to the corresponding 
extreme value attained (±1) and their dynamics stopped. These agents cannot modify their state in 
subsequent times, but they are still linked to the network and taken into account in the dynamical 
evolution of the undecided agents. In this sense such limit values can be interpreted as final states of 
irrevocable decision. 

The rewiring scheme involves cutting certain links and creating new ones, as explained in detail in fl8'. 
At the time of cutting an agent i preferentially breaks its link with agent j if there is large disagreement. 
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as quantified by pij = \xi — Xj\/2. Explicitly, its neighbours are chosen in decreasing order of the opinion 
difference Pij for p^j > 0. Then the agent creates the same number of new links based on either of two 



link-formation mechanisms known as triadic and focal closure 19 . In the former it closes a triangle 
with the "friend of a friend" if the new link can help the agent in reaching a state of total conviction (as 
measured by qij = \xi + Xj|/2), while in the latter it creates a link with a further neighbour provided 
their opinions are similar (according to Vij = \xi — Xj\/2). Like in the cutting procedure, new links are 
created in decreasing order of the opinion similarities q^j and r,jj for Qij, r^j > 0. 

In the rewiring scheme used here, agents can create links using these two mechanisms at will. In order 
to control the proportion between focal and triadic closure events we introduce a quantity y £ [0,1], 
which can be used as a stochastic decision parameter for each event. In other words, before creating 
new links an agent chooses a random number ^ € [0, 1] from a uniform distribution, and then uses the 
focal closure mechanism if ^ < y or otherwise the triadic closure mechanism. Note that focal closure 



was not considered in 18 20] for the sake of clarity, since in the absence of an external field opinion 
formation is practically dominated by close interactions with neighbouring agents. This might not be the 
case when the important feature is the external information given to all agents in the network, thus we 
will investigate the role of y here. 

In order to apply this model to our present situation we need to reinterpret the meaning of x and give 
a reasonable form for the external field representing the information given to the social group. The field 
has to break the symmetry of the model between positive and negative a;, since now the state variable 
represents the degree of the individual agent's perception of scientific knowledge. We interpret x = 1 
as expert knowledge, only attained by learned individuals in the field of science of the question or issue. 
On the other hand ignorance is represented by a very small magnitude |a;| ^ 0, whereas x < describes 
people inclined to disagree with the sound scientific information or agree with unsound concepts, e.g. 
superstition. Total opposition to scientific truth corresponds to a; = — 1 and we shall refer to such agents 
as fundamentalists from now on, since their attitude against new scientific evidence only reinforces their 
beliefs and instead of learning, they are stubbornly defending their position against. 

As the external scientific information is represented by the field hi in Eq. ([2]), we assume that hi > 
represents a drive towards scientific truth, while hi < describes a drive towards fallacy. The way this 
information is perceived by a particular agent should then be asymmetric, in the sense that it should 
help agents to learn and agree with the scientific truth. The simplest way to appropriately break the 
symmetry is by the linear expression, 

hi = h{l - Xi), (5) 

where /i is a constant strength representing the raw external information given to everybody in the 
network, and the term in brackets reflects an instantaneous personal reaction to the field. Observe that 
agents with positive x are less affected by this term, while the ignorant and fundamentalist agents are 
very much in conflict with the field, either positive or negative. This mimics the fact that superstitious 
people are in general more prepared to believe anything without proof, and change their position with 
ease. We have also tested quadratic and cubic expressions for the field and the qualitative behaviour of 
the model remains unchanged, although its analytic treatment becomes unnecessarily complicated. Thus, 
we have opted for the simpler term in Eq. ([s]). 

It should be noted that the effect of the attitude parameter is not directly affecting the response 
to the field, but it becomes extremely important when the third term in Eq. ^ compensates the second 
term, producing an exponential approach to the opinion limits dictated by the first term. In what follows 
we exhibit numerical results for this model. 

Effect of external field or "mass media" 



We have performed extensive numerical calculations using the model described above. We first adjusted 
its parameters to produce the results already presented in 18 20 , where opinions and network evolve in 
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the absence of an external field {h = 0). We initialized the system to a random network configuration of 
N = 200 nodes and average degree (/co) = 5, and kept it so for all calculations. We have tried networks 
of various sizes, and found that the effects of size scale exactly in the same way as the original model 



without field 18 . Accordingly, we fixed a set of random values for the individual attitudes ai from a 
uniform distribution with center ac = and unit half- width, and chose the initial values Xi{0) = x1 of 
the state variables from a Gaussian distribution with zero mean, unit standard deviation and cut-off at 
\x\ = 1. With this set of parameter values the network splits in communities sharing the same opinion 
for g = 10"^, where dt is approximated by a time step of size At = 10"''. 

In order to eliminate the effect of randomness in the calculations we start by setting y — Q, i.e. the 
rewiring takes place only by the triadic closure mechanism. In this way we avoid the need for making 
averages and thus probe the sole effect of the field h on the system by keeping the initial conditions fixed. 
In Fig. [l] we show the himmeli 30 visualisation of the asymptotically stationary final state of the network 
after simulating 10^ time steps of the dynamics. The right and left columns correspond to positive and 
negative field respectively and for the parameter values described above, as the magnitude of the field h 
is increased. 

These results show interesting effects of the external field on the configuration of the network. First 
of all, the ratio between the number of experts {x = 1) and fundamentalists {x = —1) grows for positive 
field and diminishes for negative field, as expected, but their distribution in the system is not symmetric. 
This means that for growing positive field quite a few fundamentalists linger in several groups, and all 
eventually join in a single community that is much more interconnected than the communities of experts. 
On the other hand, when a negative field grows in magnitude, i.e. an increasing effort to convey non- 
scientific information to the public, the expert agents opposing the fallacy are dispersed in the network 
and do not form a community. In terms of the public perception of a concept promoted by raw external 
information /i, this implies that scientifically sound concepts (associated with h > Q) require a larger 
field magnitude to create opinion consensus in the network than concepts not validated by science (i.e. 
/i < 0), since the fundamentalists found for h > Q organize themselves in close communities that prevent 
consensus of opinion. 

This surprising effect is in agreement with the behaviour of some real social networks (for example, 
creationists are well organized in very interconnected communities [31| ). Furthermore, Fig. [l] shows 
that for growing negative field the agents attain a stationary state before having time for rewiring, so 
the network remains as random as initially but with all of its agents decided already. This situation is 
asymmetric with respect to the sign of h, since community structure is preserved under a growing positive 
field, and even at high values of /i > the network's connectivity is larger than is the case for its random 
counterpart. 

Indeed, for progressively stronger positive field the number of undecided agents grows steadily, while 
for negative field there are zero agents that "do not know". These undecided agents (some of them 
fundamentalists for a weaker positive field) are forced by the large magnitude of h to have the expert 
opinion but continually resist to do so. Such undecided agents have small values (i.e. behaving as 
ignorants), and go from being distributed at the edge of the larger cluster of experts to connect the group 
of fundamentalists with the rest of the network. In this sense the ignorants act as bridges between the 
tight fundamentalist community and the well-informed people. 

In order to analyse the relationship between opinion and attitude in our model systems, we look 
at the ensemble averages and classify the agents into groups denoted by [x, a] . Here the labels take 
the symbolical values x = -1-1, -|-0, — 0, — 1 for experts, ignorants with positive opinion, ignorants with 
negative opinion, and fundamentalists, respectively, and a = — for attitude parameter > ac and 
ai < ac- In Fig. [2^ we plot the fractions /x.a of agents in each group [x, a], to show the relative group 
sizes after the normalisation condition ^ /x.q = 1. Here the thin lines are drawn to distinguish all eight 
possible groups by value of x and a, while thick lines are drawn to separate the contributions of experts, 
ignorants and fundamentalists. It is clear that a negative external field produces roughly symmetrical 
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attitude distributions in groups of decided agents, a steady asymptotic growth to negative consensus, and 
a total lack of undecided agents. On the other hand, for growing positive field it turns out that experts 
have mainly positive attitude, all fundamentalists have negative attitude, and there is a large amount of 
ignorants. 

The asymmetry due to the external field is even more evident in the actual topology of the system. 
In Fig. [2|d we show the relative contributions of the different groups to the average degree of the network 
(k) as measured by the quantity fx,akx.a/k, where ky^ ^ is the average degree of agents in group [x, a]. 
In this case we have a similar normalisation condition, fx,akx,a/k = 1. The symmetry in attitude 
sign for negative field is also seen in the left part of Fig. [2]3, but for positive field the undecided agents 
(forming the majority of the network, as shown in the right part of Fig. [2^) contribute relatively little 
to the total degree, since most connections belong to expert agents with positive attitude. Moreover, the 
fundamentalist community found for low h > has many connections when compared to its size. 

In addition to the simulations, considerable insight into the behaviour of the coupled dynamics of 
opinion and network structure can be gained with a relatively simple yet analytical mean-field approach, 
described in detail in section Materials and Methods. It turns out that Eq. ^ has an asymptotic 
solution Xi — {x^ — x*)e'^* + x* , where the fixed point x* — {h + acfi)/{h — fs) and its associated 
eigenvalue X = fs — h can be derived explicitly in terms of h and the other parameters of the model. As 
shown in Fig. [2]:, increasing the magnitude of negative field results in A > and a corresponding repulsive 
X* that grows towards a; = 1, thus giving rise to a diminishing amount of experts and building negative 
consensus. For increasing positive field the eigenvalue changes sign at a critical value ho ~ 10, signalling 
a transition to a state where an attractive fixed point hinders most agents from attaining extreme values 
of opinion. As h grows x* moves once again towards x ~ 1, implying positive consensus in the limit of 
infinite field. Fig. [2]; shows good agreement between the numerical computation of x* and A and their 
analytical approximations. 

It is also instructive to look at the effect of h on the weighted nearest-neighbours' average opinion for 
agents in group [x, a], denoted by fx^ax""^- seen in Fig. [2h, the change of sign in the experts' 
for increasing magnitude of the negative field implies that tney are sparsely distributed in the system, 
while for < h < ho the fundamentalist group is tightly interconnected and manages to keep a negative 
nearest-neighbours' average opinion until its disappearance. A phase change in the number of undecided 



agents (reminiscent of the one discussed in 32 , in that case as a function of the time-scales' ratio g) also 
shows in the value of x""^ for both groups and the network as a whole. This is signalling a very slow 
transition to positive consensus as the positive field is increased. 



Effect of focal versus triadic closure 

We now turn to explore the changes introduced by varying the parameter y from zero to one, measuring 
the relative amount of focal closure instead of triadic closure mechanisms used by agents in the network 
rewiring process. The effect of y in the final state of the opinion dynamics is illustrated in Fig. [3] for 
external field h — —10, 10. We used the same set of initial conditions and parameters as in Fig.[TJ keeping 
just the inherent randomness associated with non-zero values of y. The main difference with the y ~ 
case is due to the fact that when y grows it is easier for agents to find someone to create a new link 
with, resulting in a systematic loss of heterogeneous structure. Indeed, for /i > and as y is increased 
communities of the same opinion quickly merge, the average degree in the network and the number of 
connections between the two remaining clusters grow, and the number of ignorants decreases. In the case 
of negative field the network's structure is basically random for most values of y, except for j/ ^ 1 where 
the distribution of connections is more heterogeneous. 

There are several qualitative properties of the system that hold for all values of y. First, the number 
of experts for positive field is always smaller than the number of fundamentalists for negative field of 
the same magnitude. One could infer from this that a true scientifically sound concept is more difficult 
to acquire than a wrong concept without scientific content. Second, the agents with no neighbours 
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("outcasts", so to speak) are usually the first agents to become decided in the dynamics, either experts or 
fundamentalists, and their number grows with y. Third, the communities of fundamentalists for positive 
field are more tightly connected than the corresponding groups of experts for negative field of the same 
magnitude, or in other words, disagreeing individuals are more sparsely distributed in the network for 
h < Q than for h > Q. Finally, the undecided agents for positive field surround the largest cluster of 
experts and serve as bridges to the fundamentalist community. 

To further validate these results and take into account the relationship between opinion and attitude, 
we again perform averages over the relevant random initial conditions and parameters. In Fig.|4^ we plot 
the effect of y on the relative contributions of experts, ignorants and fundamentalists, corresponding to 
the thick lines in Fig. [2^. Here we see that the number of undecided individuals gets minimized around 
y ~ 0.7. This particular ratio between focal and triadic closure mechanisms optimizes the presence of 
experts for h > Q, although there is a persistent group of ignorants with both negative opinion and attitude 



constituting roughly 20% of the network. This is reminiscent of a result in the 2005 Eurobarometer 33 
where the validity of 13 scientific-knowledge statements (such as "It takes one month for the Earth to go 
around the Sun") was asked and the average percentage of wrong answers amounted to 21%. 

The minimum in the number of ignorant agents also shows in the relative group contributions to the 
average degree of the network, as shown in Fig. ^p. The fundamentalist group found in the < h < hg 
region grows with y, has mainly agents with negative attitude and a relatively high amount of connections 
and triangles. Quite surprisingly, the average clustering coefficient c of the network (not shown here) 
does not diminish when y ^ 0, but increases. This happens in spite of the fact that focal closure is not 
an explicit mechanism for the formation of triangles. 

Finally, in Fig. [ij; and d we can see the effect of y 7^ on the eigenvalue A associated with Eq. ^ and 
the weighted nearest-neighbours' average opinion of undecided agents /±o,±^±o±' opposed to their 
values for y = shown in Fig. [2]: and d. In Fig.|4j: we note that the critical point Hq (signalled by a change 
of sign in A) moves with due to a non-trivial increase in the eigenvalue for the region < h < ho, 
and attains a maximum value of ~ 30 for y ^ 0.7. In Fig. [4]i we observe that the position of the phase 
change in the number of undecided individuals (occurring at zero field for y = 0) is also displaced with 
increasing y, and maximized to a value close but less than Hq. Incidentally, this also happens for y ^ 0.7 
and thus coincides with the minimum in the number of ignorants after the phase change has taken place. 

In the next section we will discuss and interpret these results in the light of social behaviour, by 
adjusting the field strength h with the help of actual data extracted from extensive polls made in Mexico 
and Europe. 

Adjusting field strength with survey data 

Surveys and polls do not exactly refiect public opinion, but at least they provide some measure of the 
public perception about the subjects under investigation. The reasons for the inaccuracies are many, 
not the least how the survey was conducted including what was asked and how the questionnaire was 
designed. Most surveys could be regarded as snapshots of society without dynamical information. The 



survey carried out by the National Science Foundation 34 is an exception, where the data has been 
integrated from 1979 to 2001. Furthermore, it is difficult to follow the time development of opinion, as 

i.e. 



there has been a shift from measuring mostly literacy to matters concerning science and society 35 
from concentrating on knowledge questions to asking about attitudes and trust. 

It is quite common that the results of the survey are presented by giving percentages of the responses 
of the population to a few answering-options in the questionnaire, without further details on the topology 
of the underlying social network and of possible relations between its individuals. This lack of detail about 
social structure is unfortunate, since (as evident from above) our model could give very rich information 
about how society is organized and functioning. Irrespective of these shortcomings the surveys are 
considered useful, since we can adjust the field strength in our model to give the same percentages of a 
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population's responses as in the actual survey data, and then compare the resulting network topologies 
that would correspond to surveys of different populations. 

Here we will use two different surveys: the well-known Eurobarometer (EU) [36 1 and a Mexican survey 



(Mx) 37 . The Mexican Consejo Nacional de Ciencia y Tecnologi'a has carried out surveys on science and 
technology perception every two years, starting from 2005. The 2009 survey covered the whole country, 
with subjects aged 18 or more, giving a total of 40 469 253 citizens (of whom around 55% are women 
and 45% men). Although the Mx survey follows the methodology of the EU survey, the multiple- choice 
answers are different. In the Mx survey the possible answers to a given statement are: totally agree, tend 
to agree, tend to disagree, totally disagree and do not know, while in the EU survey there is an extra 
option: neither agree nor disagree. Thus, in order to perform a comparison between surveys we consider 
{tend to agree, totally agree} together, the same for {tend to disagree, totally disagree}, and exclude the 
do not know and neither agree nor disagree answers, since they depend too much on the protocol. The 
agreement and disagreement percentages are then normalized by its sum to give the fractions n+ and 
n_, respectively. 

The problem in starting a comparison of the model results with survey data is that there is no clear 
way for assigning an h value to each statement. The main quantitative information at our disposal is the 
percentage of agreement, and we shall use it in the following way. To assign an h value to each statement, 
we select a few statements from the EU survey and use our judgment to order them from blatant fallacies 
to obvious facts. A list of statements in this order is shown in section Materials and Methods, with 
their percentages of agreement and disagreement in Table [l] The fraction n+ is also calculated from the 
numerical results and normalized in the same fashion as with the data, which we show as a continuous 
line in Fig. [5^ and b. Since this line is mostly a monotonic function, we can give an h value to each 
statement so that the fractions n_|_ between survey and model are the same. The same procedure is then 
followed with the Mx survey, where one finds different values of the field for the same statements. The 
adjusted values of h for both surveys are shown as dots in Fig. [5^ and b. Given that our simulations of 
n+ do not depend very much on the y parameter (describing the proportion between focal and triadic 
closure), we perform the adjustment procedure by choosing only one value of y = 0.7. Notice that the 
order of the statements in the Mx survey is slightly different, and consequently the values of the field for 
a given fraction of agreement differ as well. 

Of the 15 selected statements, 9 are concerned with purely scientific matters while the other 6 also 
depend on cultural-related beliefs, although it is indeed difficult to draw a separating line between the 
two. For instance, the statement "We depend too much on science and not enough on faith" shows 38% 
agreement in EU and only 16% agreement in Mx, reflecting cultural differences, while the statement 
"Thanks to science and technology, there will be more opportunities for future generations" has 75% and 
88% agreement in EU and Mx, respectively, with practically the same proportion of disagreement. In 
Fig. [5: we compare the assigned fields of all 15 statements for the two surveys, conserving the order of 
Table T] for the EU survey and changing it for the Mx survey, as indicated by the statement labels. In 
Fig.[5p" we keep only the 9 questions without culturally-related content, shown in the same order for both 
surveys. 

First of all it is clear that total agreement or disagreement is never attained, since the dots are bounded 
by I ft I ~ 30 (an exceptional case in the Mx survey is "Smoking can cause lung cancer", where the external 
field is clearly very strong and produces 98% of agreement). A peculiar difference is that statements near 
a 50-50 agreement (i.e. ft ~ 0) are absent in the Mx survey. This we could interpret as an EU population 
less affected by propaganda modelled as a field. It is also interesting that polemic statements (for which 
the response tends to be close to 50-50) are not really distinguishable from statements lacking interest, 
for which give people give a random answer. 

Another point worth mentioning is that although our model predicts an asymmetry between positive 
and negative field, the results plotted in Fig. [5j: look more symmetric for the Mx survey than for EU. In 
terms of the description of our model, this would imply that the Mx population requires a larger external 
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influence to agree to a fallacy. Additionally, we see that for large positive field the behaviour of both 
populations is very similar, yet for moderately h > the EU survey shows a smaller response. Notice that 
for the non-cultural statements in Fig. [5]i the behaviour of both populations is also similar, from which 
we conclude that efforts to inform about matters related to science are comparable once culturally-related 
beliefs are addressed. 



Discussion 

The value of mathematically modelling a complex system, even in the case of one as intricate as society, 
is that it allows us to gain understanding of the system's response in situations not quantified in real life 
before. To this end we have modelled the effect of a controlled external field on the opinion formation 
process over a co-evolving social network, and wc have inferred its consequences in society. As the 
main feature of our model we conclude that individuals who are stubbornly opposing the spreading of 
scientific facts have to compensate it by making tightly-connected communities, where they support 
each other against commonly accepted notions. However, these communities are not completely isolated 
because ignorant agents serve as bridges connecting different parts of society. This result suggests that 
scientifically sound concepts are more difficult to acquire than concepts not validated by science, since 
opposing individuals organize themselves in close communities that prevent opinion consensus. 

We have thoroughly analysed the asymmetric response to such field in the dynamics of attaining 
complete agreement on a scientific fact or correct notion. In the spirit of a mean field approach we could 
demonstrate analytically that there is a critical field value, beyond which the dynamics is slowed down 
enormously and a heterogeneous community structure remains. This observation echoes with previous 
discussions of survey results, where it is seen that very aggressive propaganda does not necessarily result 
in a proportional immediate increase of agreement. 

We have also made an initial attempt at implementing data on scientific perception surveys from two 
different populations with our model. Even though current surveys do not fully probe the intertwined 
relationship between social communities and opinion (only measuring an averaged response emanating 
from the social dynamics), we have used their resulting agreement fractions in science-related statements 
to adjust the field strength in the model, allowing us to point out differences between the two populations. 
Particularly, we could infer that a more sceptical society has an asymmetric response in the limits of strong 
positive and negative external field, and favours divided positions for a weak external field. 

A comparison between model results and a real opinion formation process is beyond the scope of the 
present work, yet constitutes a very interesting and challenging task for the future. To advance in this 
direction one could analyse the available media content on each particular subject over time and relate 



it to a quantity like the external field h in our model, as some related studies on this matter suggest 38 
Furthermore, current surveys could be enhanced by clarifying the "do not know" answers, integrating 
data over time, and keeping track of the network topology details. 

The present study was aimed at investigating the susceptibility of an ethnic group facing such universal 
notions as science and technology. This problem has occupied the thoughts of many scholars for a long 
time, and we think it is time to make efforts in quantifying such notions. It is expected that there are 
cultural differences in answering specific questions, but it is also not surprising to find many similarities 
due to human cognitive processes and the globalized world we live in. Further insights coming from 
mathematical modelling, such as those presented here, could be useful in improving current policies to 
communicate the facts and findings of science in a socially responsible way. 
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Materials and Methods 
Analytical calculations 

We can introduce a series of analytical approximations to shed light over the complex dynamics found 
in our co-evolving social network model, and help in the understanding of the numerical simulations 
discussed in the Results section. In order to accomplish this, we analyze the stationary state resulting from 
the coupling between opinion dynamics and network structure in the spirit of a mean-field approximation. 
Let us start by introducing a site-dependent integrating factor, 

lit) = e-/'(''-/=)''^ 

for the equation of motion dtXi + (h — fs)xi — aifi + h, which allows us to write the formal solution, 



f I{T){a,fl{T) + h)dT 







(6) 



where 1° = 2^(0)- In the asymptotically stationary final state of the system, we can substitute the site- 
dependent short- and long-range interaction terms in Eqs. ( [4| and ([s]) by average quantities fs and /; 
that vary slowly in time, so that I — exp [{h — fs)t\ and Eq. ^6]) can readily be integrated to give, 

= (xO-x*)e^*+x*, (7) 

where x* is the mean fixed point of Eq. ([2| and A its associated eigenvalue. 
Indeed, a linear stability analysis of Eq. ([2| around x* results in, 

X* - x-\ - ^ f8) 



and, 

A 



d{dtXi) 



dxi 



fs - K (9) 



where the quenched attitude parameter has been substituted by its average ac and correlations were 
ignored. These expressions allow us to determine the qualitative behaviour of the system as a function of 
the external field. Since fs and /; are respectively bounded by the degree k and by iV — fc — 1, for ^ 
the eigenvalue is trivially positive, the repulsive fixed point approaches x = 1 and a steady asymptotic 
growth towards negative consensus is achieved. For h ^ this situation is reversed to give a negative 
eigenvalue, an attractive fixed point approaching x = 1, and positive consensus in the limit h ^ oo. Since 
fg is continuous, there is at least one value /iq such that /s(/io) = fiQ, an implicit definition of the critical 
point discussed in the Results section. Note that consensus is eventually achieved in both directions, but 
a fundamental asymmetry arises in the dynamical evolution of the system due to the nature of its fixed 
point, with only decided agents for negative field and undecided agents for positive field. 

The idea now is to find appropriate expressions for fs and // in terms of x* and the parameters of 
the model, then solve for the fixed point and obtain the explicit dependence on the field h. To simplify 
further derivations it is useful to introduce the scaled fixed point u — box*, where bo — •\/2/7r/erf(l/%/2), 
and with which the Gaussian distribution of initial opinions can be written as Pq{x) = {bo/2) exp(— a;^/2). 

As an intermediate step we estimate the change in the average degree of the network due to the 
first parallel rewiring of undecided agents. As described with detail in |18|, the simultaneous rewiring of 
neighbouring nodes can lead to the net creation or deletion of a link. For simplicity we assume y = (i.e. 
only triadic closure events) and fco <C N/2, so that the network to be rewired is essentially a Cay ley tree 
with coordination number ko. Since the number of 2nd neighbours of a rewiring agent i is n2 = ko{ko — 1) 
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and thus larger than its degree, agent i chooses to cut all of its fco connections and create /cq links with 
2nd neighbours of appropriate opinion. Each rewiring neighbour takes the same decision with probability 
1, thus bonds are cut twice leading to an increase in the degree proportional to fep. On the other hand, 
each newly-rewired neighbour can also choose agent i from its own set of 2nd neighbours to create a 
link with probability fcg /ri2 , thus bonds are made twice leading to a decrease in the degree proportional 
to fcg/n2- This situation requires the simultaneous rewiring of two undecided agents, then the resulting 
change in degree is estimated as fc — /cq = ?iQfco[l — l/(fco ~ 1)]: with uq the fraction of undecided agents 
in the network. 

Now, for h < and in the presence of a repulsive x* , most agents get decided before the previously 
discussed first rewiring takes place, so the calculated degree is a good approximation for the final degree 
of the network. As for no, it can be estimated by integrating the distribution of initial opinions po over 
the interval [Xg , Xg ], where xf=x*±{lT x*) exp{-\gAt) are the positive/negative initial opinions 
required to reach extreme opinions in exactly g time steps of size At, according to the solution in Eq. ([t]). 
The resulting integral can be expanded to first order for x^ < and gives no = bo exp(— AgAt). Thus, 
considering that /s ~ for early times we can write the final degree in the system as. 



1 



^ ko-ll ° 



(10) 



implying an exponential decrease in fc — /cq &s the magnitude of ft, < increases. 

Let us continue and find asymptotically correct expressions for fg and /; in terms of the scaled fixed 
point u = box* and the parameters of the model. We denote by {si Sf} the fraction of agents that start 
with an initial opinion sign Si and end up in an extreme value of opinion with sign s/ . In the presence of the 
repulsive fixed point x* found for h < 0, the final fraction of experts is {++} = po{T)dT = (1 — u)/2, 

and {H — } = po{T)dT = u/2 is the fraction of initially positive agents that convert to fundamentalists 
due to their adverse neighbourhood. On the other hand, { — h} = and the fraction of initially negative 

agents with a favouring neighbourhood is { } = 1/2. Then, the short-range interaction term can be 

written as fs = {++}f^ + ({ } + {H — }) f^^ where the corresponding contributions by experts and 

fundamentalists are, 

if+ =2k[{++}-{+-}] = k{l-2u) 
\f- =2k[{— }-{-+}] = k 

that is, 

f^ = k[u^ -u + 1]. (11) 

Here we refiect the fact that, in average, the field-favoured negative agents end up as fundamentalists 
being connected mostly between themselves, while most positive agents are in conflict with the field and 
cannot find proper neighbourhoods. 

Wc then follow a similar argument for the long-range interaction parameter and write // — {H — h}/;^ + 

({ } + {H — }) /; , where positive agents contribute with /j+ = 0^.(1 — 2m) and negative agents with 

fi^ = — Ofc, since // does not have a sgn(a;i) term like fs and consequently it is asymmetric with respect 
to opinion sign. The factor Ofc measures the average number of further neighbours weighted by their 
geodesic distance to an arbitrary node i, in accordance with Eq. ([3|, and can be easily calculated for an 
acyclic network as — k{k— lY~^ /£ + R/£maxi where £rnax — 1 is the largest layer number such 

that the sum over i is less than N, and R is the remaining number of nodes in the network, excluding 
agent i. Thus, /; is estimated as, 

/, = -ak [2u - u^] . (12) 



Finally, we can use Eq. (10) to write k and at explicitly in Eqs. (11) and (12), then insert these in 
Eq. ^ and solve for u. After some algebra we find a cubic equation for the scaled fixed point, 



ku^ + (fogckcifc — k)u^ + {k — h ~ 2boacO'k)u + boh ~ 0, 



(13) 
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which can be solved analyticahy to give x* = u/6o as a function of h and the parameters of the model. The 
result, plotted as a continuous line in the bottom part of Fig. [2]:, is in close agreement with the numerical 
simulations. With this solution we can also calculate analytically the eigenvalue X = k [u^ — u + l] — h, 
shown in the top of Fig. [2j:, and the mean value of the nearest-neighbours' average opinion x"" — 

[{++}/+ — ({ } + {H — }) /7] /k — V? — 2u, depicted as a continuous line in Fig. [2]i, again in good 

agreement with numerical results. 



Survey data 



We chose 15 equivalent statements from the EU 36 and Mx 37 surveys concerning the public perception 
of science and technology in society. Statements in both surveys are deemed equivalent apart from minor 
changes in wording and differences in country names, and according to [36| they read as follows: 

a Thanks to scientific and technological advances, the Earth's natural resources will be inexhaustible. 

b Science and technology can sort out any problem. 

c Science and technology cannot really play a role in improving the environment, 
d Science should have no limits to what it is able to investigate, 
e We depend too much on science and not enough on faith, 
f Some number are especially lucky for some people. 

g Scientists should be allowed to experiment on animals like dogs and monkeys if this can help sort 
out human health problems. 

h The benefits of science are greater than any harmful effects it may have. 

i Because of their knowledge, scientists have a power that makes them dangerous. 

j Science makes our ways of life change too fast. 

k The application of science and new technologies will make people's work more interesting. 

1 Compared with research carried out and funded by each Member State, to what extent do you think 
that collaborative research across Europe and funded by the European Union is in the national 
interest? 

m Even if it brings no immediate benefits, scientific research which adds to knowledge should be 
supported by Government. 

n Thanks to science and technology, there will be more opportunities for future generations. 

o A scientific discovery is in itself neither "good" nor "bad" , it is only the way the discovery is used 
which matters. 

According to our subjective judgment, the previous list is ordered from blatant fallacies to obvious 
facts. In Table [T] we show all statements for both the EU and Mx surveys. The percentages of agreement 
(yes) and disagreement (no) are also shown, such that the remaining percentage corresponds to the sum 
of the do not know and neither agree nor disagree answers for the EU survey, and to the do not know 
answer for the Mx survey. 
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Figure 1. Final state of the network for y = 0. Asymptotically stationary state of the dynamics 
for chosen initial conditions and parameters, as described in the text. Decided agents are represented by 
red {xi — 1) or blue {xi — —1) circles, and undecided agents by yellow (0 < a;^ < 1) or black 
(— 1 < a^i < 0) squares. The right (left) column corresponds to positive (negative) field h of increasing 
magnitude. All calculations were done with y = 0, i.e. rewiring with triadic closure mechanism only. 
The visualisation shows the basic asymmetry of the model: for stronger negative field the community 
structure is quickly lost and there is an asymptotic growth towards negative consensus, while for 
growing positive field fundamentalists linger in well-connected communities and ignorants slow down 
the drive towards positive consensus. 
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Figure 2. Effect of external field, (a) Relative group size as a function of the external field h, 
where each group [x, a] is distinguished in terms of both opinion and attitude, as described in the text. 
Thin lines separate groups and thick lines divide the contributions of experts, ignorants and 
fundamentalists, (b) Relative group contribution to the average degree of the network as a function of 
h. (c) Fixed point of Eq. ^ (bottom) and associated eigenvalue (top) as a function of the external 
field. Numerical results are shown as dots, while the corresponding analytical approximations are 
depicted as lines, (d) Weighted nearest-neighbours' average opinion for agents in group [x, a] as a 
function of h, shown as symbols for different opinion groups. The continuous line depicts the analytical 
approximation of its mean value over all agents. All numerical calculations are averaged over 108 
realisations of the dynamics with y = 0. 
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Figure 3. Final state of the network for y 7^ 0. Asymptotically stationary state of the dynamics 
for the same initial conditions and parameters as in Fig. [l] Decided agents are represented by red 
{xi = 1) or blue {xi = —1) circles, and undecided agents by yellow (0 < Xi < 1) or black (—1 < < 0) 
squares. The left and right columns correspond to h = —10, 10, respectively, as the parameter y is 
varied from zero to one. An increasing amount of focal closure events results in a systematic loss of 
heterogeneous structure in the network, higher degree, and decreasing amount of ignorants. 
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Figure 4. Focal versus triadic closure. Relative group size for experts, ignorants and 
fundamentalists (a), relative group contribution to the average degree of the network (b), eigenvalue 
associated with Eq. ^ (c), and weighted nearest- neighbours' average opinion for undecided agents (d), 
all plotted as functions of the external field h and for rewiring parameter values y = 0.1, 0.7, 1. All 
calculations are averaged over 108 realisations of the dynamics. Non-zero values of the rewiring 
parameter retain the qualitative picture of the y = case shown in Fig. [2] along with a displacement of 
the critical point ho and non-trivial behaviour in the region < h < Hq. 
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Figure 5. Adjusting field strength with survey data, (a) External field strength ft, as a function 
of the fraction of agreeing individuals n+ for a rewiring parameter value oi y = 0.7. Continuous lines 
show the model simulations for the same set of parameters as used in Fig. |4] Symbols correspond to n+ 
for statements in the EU survey, with adjusted h values so that the agreement fractions between 
statement and simulation are the same, (b) The same but with symbols corresponding to statements in 
the Mx survey, (c) Assigned h values for all 15 statements in the EU and Mx surveys, each set in 
different order according to the statement labels, (d) The same but only for 9 selected statements 
without cultural content, shown in the same order for both surveys. 
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Table 1. Survey data 





EU survey 


Mx survey 


statement 


yes (%) 


no (%) 


yes (%) 


no (%) 


a 


21 


56 


33 


63 


b 


22 


57 


21 


76 


c 


24 


54 


15 


81 


d 


35 


44 


70 


17 


8 


38 


34 


16 


84 


f 


40 


35 


34 


64 


g 


44 


37 


57 


41 


h 


46 


20 


69 


23 


i 


53 


24 


57 


39 


j 


58 


22 


83 


14 


k 


61 


14 


81 


15 


1 


66 


10 


94 


3 


m 


72 


9 


88 


9 


n 


75 


8 


88 


9 





78 


7 


91 


6 



Statements from the EU and Mx surveys that have been considered for this study. The statement labels 
are ordered from blatant fallacies to obvious facts, according to our judgment. The percentages of 
agreement (yes) and disagreement (no) for each statement and survey are also shown. 



